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We investigate instabilities in the charge channel in the vicinity of (meta-)magnetic transitions 
of itinerant electron systems. Based on a weak coupling analysis we argue that in a one-band t-t 1 
Hubbard model near the van Hove filling and dominant ferromagnetic fluctuations it is difficult to 
account for a microscopic mechanism for a d-wave Pomeranchuk deformation of the Fermi surface. 
A similar deformation has been considered for the metamagnetic transition in Sr3Ru207. As an 
alternative we discuss the possibility of charge inhomogeneity on the nano scale. This extends the 
analogy of the metamagnetic transition to a liquid-gas transition. 
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Over the last decade, the perovskite strontium- 
ruthenates have become a model system for the study of 
correlated electrons. While the infinite-layer compound 
SrRu03 is an itinerant ferromagnetQ with Curie temper- 
ature of Tc = 165K, the single-layer Sr2Ru04 is a triplet 
superconductor^ below T c = 1.5K. Above T c it is a well- 
studied strongly anisotropic Fermi liquid j3J]. By substi- 
tuting Sr with La and thereby enhancing the Fermi sur- 
face (FS) volumes, it can be brought into the vicinity of 
ferromagnetism (FM) 4]. The double-layered Sr3Ru20y 
interpolates between these two limiting cases in an ex- 
citing way. While first experiments found FM, subse- 
quent research|5j made clear that the magnetization M 
is zero in the absence of a magnetic field B. However at 
a direction-dependent critical field strength B m between 
4 and 8T, M(B) shows a steep increase below T ~ 2K, 
referred to as a metamagnetic transition 0, Q- Around 
B m , there are indications for non- Fermi liquid behav- 
ior, e.g. in the specific heat and the T-dependence of the 
inplane-resistivity0,0|. Theoretically, the metamagnetic 
transition with a single jump in the magnetization can 
be understood in a Stoner-type model with a FS near a 
van Hove (VH) singularity in the density of states 0,0. 
Indeed a number of experiments Q for the single-layer 
system S^RuCU and band structure calculations f° r 
both Sr2RuC>4 and Sr3Ru2 0y confirm the presence of a 
nearby van Hove singularity for the FS derived from the 
Ru d Xi) -orbitals. 

A new twist occurred when improved sample quality 
revealed at least two iumns|llL Il2| in the magnetization 
when B is increased. Most interestingly, in the critical 
field range between the two jumps of Ref. 0j, the resid- 
ual resistivity po is about twice as high as below and 
above the jumps. As a function of the temperature, p(T) 
saturates already at T = IK in the critical field range, 
and the T-dependence above is of non-Fermi-liquid type. 
The high po can be understood if one assumes the for- 
mation of domain- like structures in the critical range [l3j. 
Recently, two works [3 El suggested the occurrence of 
a d,; 2 _ y 2 -wave-like Fermi surface deformation (dFSD). In 
a one-band picture on the two-dimensional (2D) square 



lattice, the <iFSD splits the two inequivalent VH points 
and thereby gains energy if the forward scattering in the 
charge channel has an attractive d-wave component. This 
possibility has been analyzed for Hubbard-like models 
by a number of authors |l5l llfil Il7l Il8l Il9| . An experi- 
mental realization of the effect has not been found yet. 
For the double-layer ruthenate this picture is however 
quite appealing 13]. The occurrence of two magnetiza- 
tion jumps can be understood in the dFSD framework 
as first-order entry into and exit out of the phase with 
reduced FS symmetry (lif. Moreover the sensitivity of 
the experimental picture with respect to sample quality 
seems to fit into the scenario of an unconventional order 
parameter. We note that while the chcmical-potential- 
tuned quantum phase transition into the dFSD state is 
typically first order, the thermal transition at T c is second 
order over a wide parameter range in the simple forward 
scattering model |2C)l l2l| . 

So far the microscopic origin of an attractive ci-wave 
component in the forward scattering has not been ad- 
dressed in the context of the metamagnetic transition. 
In this work we analyze this question based on the one- 
band t-t' Hubbard model on the 2D square lattice and 
specify the parameter region in which the dFSD is favor- 
able. Furthermore we reinvestigate the Stoner model for 
the metamagnetic transition 0. We point out that in this 
model, the magnetization jump is located in an unstable 
density region. This extends the analogy to a liquid- 
gas transition already emphasized by Green et al.[22j. 
Around the transition charge inhomogeneity might oc- 
cur and give rise to a higher residual resistivity. We also 
describe a possibility how more than one jump in the 
magnetization can occur in this model. 

First we consider the effective forward scattering in the 
charge channel for the 2D t-t' Hubbard modeljl^ on the 
square lattice in zero magnetic field. A perturbative but 
unbiased way to arrive at effective low-energy or low-T in- 
teractions is the functional renormalization group (fRG). 
In fact, the <iFSD or Pomeranchuk tendencies in the 2D 
Hubbard model were revealed using fROjlfij. Here, we 
use the temperature- flow fR.c[24l|. It has the additional 
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advantage of not being blind with respect to long-wave 
length particle-hole fluctuations. Using this method, a 
low-T p-wave superconducting instability was found for 
a FS similar to that of the 7-band in S^RuC^j^. A 
second order transition into a state with broken sym- 
metry is indicated by a divergence of the corresponding 
susceptibility at T c when T is lowered. 

The fRG flow is started at high temperatures T ~ 
bandwidth with pure onsite repulsion U . In the flow 
to lower T, one-loop corrections generate a detailed de- 
pendence of the interactions Vr(fci, ki, ka) on the wave 
vectors fci, hi 1 and k%. In the Hubbard model with 
U = 3* near the VH filling, the fRG basically finds 
two regimes 24] . One is for \t'\ < *|/3, where the lead- 
ing ordering tendencies are either in the antiferromag- 
netic channel for strong nesting, or, if the nesting is 
weaker, in the d x 2_ y 2-wa,ve pairing channel |24|. For for 
|*|/3 < |*'| < 1/2, the leading instability is in the FM 
channel. With bare onsite interactions the charge chan- 
nel never becomes leading. However, e.g., the forward 
scattering V r (k,k')cr cV, ,cr, ,cr with zero momen- 

L \ ' / /c,s /c ,s k ,s k,s 

turn transfer can develop a pronounced dependence on 
k and k' . This can give rise to subdominant instabilities 
such as FS deformations The suggested d-wave FS 
deformation needs a component in V c (fc, k') of the type 
—V c> 2(cosk x — cosky) ■ (cos £4 — cosk' y ) with V c ,2 > 0. 

In Fig. H a ) w e show the fRG result for V c (k, k') at low 
temperatures above the runaway flow of the leading ten- 
dency, obtained with the fRG scheme of Ref. |3 The 
dashed line is for temperatures above a d-wave pairing 
instability at the VH filling for t' = 0.25* and T = 0.09*. 
One can clearly observe the attractive d-wave component 
in V c (k, k') with k' fixed at (— tt, 0). The solid line is near 
VH filling and t' = 0.44* in the FM regime at T = 0.04*. 
Now there is no sign for an attraction in the d-wave chan- 
nel. The couplings to external static and uniform sources 
in this channel does not increase for T — ► (Fig. Hb)). 
These results reflect the influence of the spin fluctuations 
on the forward scattering. Quite generally, the forward 
scattering between the parts of the FS connected by the 
dominant spin fluctuation wavevector q is renormalized 
to stronger repulsion. The strong (tt, 7r)-fluctuations be- 
tween the VH points create the +- lobe of the d-wave in 
V c (k, k') in the AF/d-wave pairing regime for small *'. 
The attraction for k sa k' is a precursor of the pairing 
tendencies. In the FM regime, the dominant spin fluc- 
tuations have q « 0, i.e. make k ss kl more repulsive. 
Hence the d-wave charge forward scattering turns out to 
be slightly repulsive rather than attractive. In fact, there 
is not any clearly attractive component in V c (k, k') in the 
FM regime. When we move the FS further away from the 
VH points by a density change or by a weak Zeeman split- 
ting, the FM instability is cut off. Nevertheless, the main 
features of the charge forward scattering remain robust. 

We conclude that at least in a single-band scenario 




FIG. 1: (color online), a): Effective forward scattering in the 
charge channel V c (k, k 1 ) with k! fixed near (— 7r,0) (6 — —tt) 
obtained with 96-patch T-flow fRG. The inset shows the Fermi 
surfaces, b) : T-flow of the coupling to static external dFSD 
fields, averaged around the FS. In all plots, the solid lines are 
for t' = OA-it, the dashed lines for t' = 0.25t. 

with dominant FM fluctuations, it is difficult to iden- 
tify a microscopic mechanism for a d-wave FS deforma- 
tion. Conversely, if we assume that the relevant FS of 
Sr 3 Ru 2 07 is on the AF/d-wave side, the spin response 
should be strongest near the wavevector q = (tt, tt) con- 
necting the two VH regions. This seems unlikely in view 
of the strong FM tendencies in all layered strontium- 
ruthenates, and in inelastic neutron scattering there is 
no sign for these fluctuations 25]. Furthermore, at least 
in the Hubbard model, the dFSD tendencies are always 
weaker than other instabilities. Hence, we would need a 
mechanism that suppresses the leading tendencies. 

We now turn to a second scenario for the metamag- 
netic transition that exploits the vicinity to ferromag- 
netism. This tendency is common all the strontium- 
ruthcnatcs mentioned in the introduction. In contrast 
with the dFSD instability it fits readily into the pic- 
ture drawn by a fRG treatment, as FM is found as the 
leading instability of the weakly-coupled Hubbard model 
near the VH filling for 1/3 < |*'| < 1/2 24]. For this 
situation, the metamagnetic transition was modeled by 
Binz and Sigrist 8] using a mean-field theory with a den- 
sity of states p(e) = ^y\og{W/\e — 6vh|] in a band 
—W<e~ evH < W. The electrons interact with an 
onsite interaction U > 0. Their density n and the dis- 
tance to the VH filling at n = 1 can be adjusted by eyH- 
Minimizing the Gibbs potential G(T, h, n) with respect 
to the magnetization m yields a density region near the 
VH filling at low T where m exhibits a sudden jump at 
a critical magnetic field h m . Coming down from high T 
at h m , to increases strongly and finally saturates below a 
temperature scale T m . Impurities cut off the logarithmic 
growth of the DOS near the VH filling. Hence, they are 
expected to smear out the metamagnetic transition. In 
this simple model, the magnetization does not depend on 
the direction of the /i-field, but this effect can be taken 
into account, e.g., by including anisotropic interactions. 

In Fig. 0a) we show the steplike increase of the magne- 
tization to with increasing magnetic field h for different 
electron densities n, equivalent to the data of Ref. In 



3 



Fig. b) we show the derivative of the Gibbs potential 
G(T, h, n) with respect to n versus n for different h. For 
the critical h- values h rn where the jump in m occurs, the 
curvature of G(T, h, n) with respect to n, corresponding 
to the inverse compressibility, is negative. Hence, at the 
metamagnetic transition, the system is unstable with re- 
spect to phase separation. If we plot G(T, h, n) versus n 
we find a small upward dent in the critical density range. 
In this range we can minimize the energy by the Maxwell 
construction, i.e., by bridging the dent from below with 
a straight line which lies below G(T, h, n) between n < (h) 
and n > (h). The average density is a mixture between a 
less dense phase n < (h) with lower m and a denser phase 
n>(/i) with higher m, the mixing ratio < p < 1 is 
determined by pn < (h) + (1 — p)n>(h) — n. When we 
increase h from below at a fixed total density n, we enter 
the unstable region at a certain lower critical field h < (n) 
with n < (h < (n)) = n. In Fig. HJb) for n = 0.94 this oc- 
curs for h a little higher than h — 0.0005TT for curve 
2. For curve 3 the stable solution is a mixture between 
n < (h) and n > (h). The path of n < (h) and n > (h) with 
increasing h is shown by the dashed lines. Going up, the 
weight p of n < (h) decreases until it vanishes for curve 4. 
In the coexistence region, fi = dG/dn is constant and the 
inverse compressibility k _1 oc d 2 G/dn 2 vanishes. With 
increasing h, the instability region moves to smaller den- 
sities, and above an upper critical field h > (n), the ho- 
mogeneous density n — 0.94 has the lowest Gibbs energy 
again (curve 5). The wiggle in dG/dn, i.e. the energy 
difference between the phase separated and the homoge- 
neous solution, becomes smaller with increasing h. 

Again, this is very similar to a liquid-gas transition. 
There the intensive tuning parameter is the temperature 
T instead of h. At the transition, the conjugated exten- 
sive quantity, the entropy jumps upwards. As a function 
of the volume V, the Gibbs potential has negative cur- 
vature around the critical volume where the transition 
occurs, and the isothermal lines in the pressure (P = 
— dG / dV)-Yo\\xme diagram exhibit a similar wiggle as 
fi = dG/dn in our case. Hence the inhomogeneous phase 
at the metamagnetic transition just corresponds to the 
coexistence regime in the liquid-gas transition. 

As can be seen Fig. b), in our model the density 
difference between the two components of the mixture 
is rather weak of the order 0.01/site. In a more real- 
istic model including long range Coulomb interactions, 
the phase se par ation will be frustrated on a microscopic 
length scale 26] . Then the two-phase mixture will cre- 
ate a pattern of nano-domains. Elongated domains like 
stripes are a likely candidates. Here we speculate that 
this charge inhomogeneity may be responsible for the 
observed plateau in the residual resistivity at the meta- 
magnetic transition. The multiple Fermi surfaces of 
Sr 3 Ru207 might help to reduce the Coulomb energy of 
the charge distribution by mutual screening. 

In the liquid-gas transition the entropy jump is con- 



a) b) 




FIG. 2: (color online), a): Magnetization m at tempera- 
ture T=0 vs. magnetic field h in units of the half bandwidth 
W and density n, U — 0AW. The two thick lines bound 
the phase separated region of the Maxwell construction, b): 
-H = -dG/dn for h = 0.0002 . . . 0.0019VK with h increasing 
from curve 1 to 5. The dotted lines indicate the Maxwell con- 
struction and span the phase separated region in the upper 
left plot. The dashed lines show the density components vs. 
h for average density n = 0.94. c): m at T=0 vs. h for total 
density n = 0.94. The thin line is for the homogeneous case, 
the dashed line is the Maxwell construction allowing for zero 
interface energy Gi = 0. The thicker line with two jumps is 
for Gi = 5 • W~ 7 W. d): Boundaries of the inhomogeneity re- 
gions for total density n = 0.94 vs. h and T. The dotted lined 
is again for Gj = and the thick line for Gi = 5 • 10~ r W. 

served by the Maxwell construction if the transition is 
crossed at fixed pressure. In the metamagnetic case the 
total density is fixed during the transition. The Maxwell 
construction interpolates between the low magnetization 
below the transition and the high-m phase above. This 
transforms the step at h m into a continuous rise be- 
tween the lower and upper coexistence fields, /i< and h > . 
Across the transition, a fraction < p < 1 will be already 
in the dense large-ft. phase corresponding to n > (h), while 
1 — p will remain in the \ow-h less dense phase with n < (h), 
and p increases from to 1. Now, the Coulomb interac- 
tion will enforce charge neutrality on a nano scale. It 
is natural to assume that the interface between the two 
phases and the remaining Coulomb energy due to the in- 
homogeneity costs some energy Gj > 0. Taking this into 
account, the straight line for the energy of the phase- 
separated solution which undercuts the Gibbs potential 
of the homogeneous system between n < (h) and n>(h) 
gets shifted upward by Gj. Then nano-scale phase sepa- 
ration is energetically favorable only for p- values between 
Pmin and pmax, otherwise the homogeneous state will be 
lower in energy. When we compute the magnetization 
with this additional constraint, p m in < p < p m ax, we ob- 
tain two steps. One is near h <1 where p jumps from 
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to a nonzero p min . The other is near h > , where the mi- 
nority domains are squeezed out, and p jumps from p max 
to 1. Hysteresis is most likely at these two jumps, in 
between the evolution is a continuous change of the rela- 
tive size p of the two components and their densities. A 
precise calculation of Gi is difficult. In order to outline 
our idea, in Fig. Elc) we assume a value G[ = 5 • 10~ 7 W 
which is about 1/5 of the maximal energy gain of the 
inhomogeneous solution. Obviously, if Gi is too large, 
inhomogeneity does not occur. 

In the experiments^!, the resistivity anomaly be- 
comes smaller when the in-plane component of the mag- 
netic field is increased. Just as was suggested for the 
dFSD domains 0, this could have the effect of orient- 
ing the density pattern. Then the transport along and 
orthogonal to the in-plane field should be different. 

When the temperature is increased the magnetization 
jumps get smoothed out. As the spin symmetry is bro- 
ken by h ^ 0, there can be no phase transition associ- 
ated with a symmetry breaking. However the instability 
towards phase separation described above yields a ther- 
modynamical distinction between the low-T micro-phase 
separated state and the normal state at higher T. Above 
a threshold temperature T neg , the negative curvature re- 
gion in the Gibbs potential disappears. Hence the density 
remains homogeneous through the /i-range where m rises 
steeply. We speculate that T neg is an upper bound for the 
temperature scale below which the T-dependent inelastic 
scattering becomes unobservable in the experiments. If 
the interface energy Gi is T- independent, the inhomo- 
geneity only sets in at T < T ncg . In Fig. Eld) we show 
the inhomogeneity regions in the h-T diagram for Gj = 
and for a T-independent Gi = 5 • 1Q~ 7 W. Lifetime ef- 
fects due to impurities act like increasing T and reduce 
the inhomogeneity region as well. 

Experimentally, this scenario could be checked by local 
probes which are able to resolve spatial variations of the 
charge distribution, e.g. scanning tunneling microscopy. 
Note that in Ref. 13] the authors state carefully that 
phase separation cannot be ruled out near the transition. 

In conclusion, we analyzed two possible scenarios for 
the metamagnetic transition in the layered ruthenates. 
For the d-wave FS deformation picture [l3l fl4|. a simple 
one-band t-t' Hubbard model can account for a charge 
instability that splits the saddle points of the dispersion 
only if the dominant magnetic fluctuations have a large 
wave vector connecting these saddle points. This does 
not seem to be the case in view of the many experimen- 
tal and theoretical indications that strontium-ruthenates 
are close to ferromagnetism. For dominant FM fluctua- 
tions, we do not find an attractive coupling constant for 
the cZ-wave FS deformations. As an alternative we ex- 
tended the Stoner model of itinerant metamagnetism by 
Binz and Sigrist and analyzed the possibility of micro- 
scopic phase separation. The metamagnetic transition 
in this model occurs in a region which is unstable with 



respect to density variations. Due to this general ten- 
dency at first order transitions the charge distribution 
might become inhomogeneous. We discussed the effects 
of the Coulomb repulsion and interface energies. If they 
do not cancel the energy gain by the phase separation 
completely, nano-scale inhomogeneity could account for 
the observed anomaly in the residual resistivity and the 
occurrence of two steps in the magnetization. 

Finally we note that Binz et al.jl^l have considered a 
scenario with uncharged magnetic domains due to the de- 
magnetization effect. We hope that further experiments 
can distinguish between the various proposals. 
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